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C<l . Abstract 

i. 

The molecule-like structure of the C isotopes (A=12, 14, 16) is investigated 
OO \ using a microscopic a+Q:+a+n+n+- • • model. The valence neutrons are 

classified based on the molecular-orbit (MO) model, and both yr-orbit and 
(j-orbit are introduced around three a-clusters. The valence neutrons which 
1^ . occupy the vr-orbit increase the binding energy and stabilize the linear-chain 

\ of ?>a against the breathing-like break-up. However, ^^C with the vr-orbit 

\ does not show clear energy minimum against the bending-like path. The 

■ combination of the valence neutrons in the tt- and the cr-orbit is promising to 

stabilize the linear-chain state against the breathing- and bending- modes, and 
it is found that the excited states of -'^^C with the (3/2~)^(l/2~)^ configuration 
for the four valence neutrons is one of the most promising candidates for such 
structure. 
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I. INTRODUCTION 

A survey of the molecule-like structure is one of the most challenging subjects in light 
neutron-rich nuclei. In Be isotopes, recently, decay into fragments of He isotopes (^He, ^He, 
^He) has been observed from the excited states of ^^Be [|1| and ^^Be [|2|,|^, and the presence of 
a two-center configuration is suggested. From the theoretical side, these states are studied by 
various models by which the molecular-orbital nature of the weakly bound neutrons around 
the two a-clusters has been revealed 0-^]. 

As for multi-cluster configurations in light a-nuclei beyond the two-center systems, the 
existence of Na states has been predicted around the threshold energy in the so-called Ikeda 
diagram P]. For example, it has been suggested that the second 0+ state of ^^C has a 3a-like 



molecular configuration [|T0|. However, according to many theoretical analyses [[T^], the state 



is not necessary to have a linear-chain of 3a, but is described as a weak-coupling state with 
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a triangular shape or the ^Be+a configuration. Furthermore, a rotational band of ^^O with 
a very large moment of inertia has been observed around the 4a threshold energy region 



T2[, and the 4a linear-chain have been discussed [|I3|,|I4[. However, the 4a-chain has not 
been experimentally confirmed yet. 

Recently, the discussions of the well-developed cluster structure are extended to the 
neutron-rich nuclei, and the role of valence neutrons which stabilize the linear-chain structure 
has been pointed out. For example, von Oertzen has extended his analyses for the molecular 
structure in Be isotopes fl^ to C isotopes, and the linear-chain state consisting of 3a and 
valence neurons around it has been speculated. Even if the 3a-system without valence 
neutrons (^^C) does not have a linear-chain structure, the valence neutrons around it are 
expected to increase the binding energy and stabilize the linear-chain state. 

In this paper, the molecular-orbital approach introduced for the Be isotopes is applied to 
a study of C isotopes (A=12, 14, 16). The appearance of prolonged structure with the a-a 
core in -^^Be and -^^Be [|1^] has been successfully described by the a+a+n+n+-- model, 
where the orbits for the valence neutrons are classified based on the molecular-orbit (MO) 



model |T^. In the present study, MO is extended to C isotopes, we focus on the stabilization 



of the linear-chain structure with the 3a-core to which valence neutrons are added. 



II. MODEL 



The total wave function of a microscopic a+a+a+n+n--- model is fully antisymmetrized 
and expressed by a superposition of basis states centered to different relative distances 
between the a-clusters (rf) with various configurations of the valence neutrons (cl, c2 ••) 
around the a-clusters: 

^iK= E ^i^i.^[^S"vrvrH0fxi)(02'x2)--]. (1) 

d,cl,c2-- 

The projection to the eigen-states of angular momentum (Pmx) performed numerically. 
Each a cluster consisting of four nucleons is described by Gaussians (Gai) centered at R^i 
and spin-isospin wave function(x): 

^(a) ^ G^^lptp'^lpllxp.XplXn.Xni, ^ = 1, 2, 3. (2) 
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l/2/?^ 



(3) 



where, the oscillator parameter {(3) is equal to 1.46 fm. For the linear-chain state, the values 
of {Rai} are — 0, and +d on the z-axis. Each valence neutron (0f Xi) around the a-a-a 
core is expressed by a linear combination of local Gaussians: 



(4) 



These valence-neutron orbits are classified according to the MO picture [|T^]. The orbit of 
the valence neutron perpendicular to the z-axis of the 3a linear-chain core is called vr-orbit. 
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and one along the z-axis is called a-orbit. The antisymmetrization imposes the forbidden 
space for the valence neutrons; the 7r-orbit must have at least one node perpendicular to the 
z-axis, and the a-orbit must have at least three nodes since three a-clusters along the z-axis 
already occupy the orbit als with = 0, 1, 2. The amplitude of the tt- and the a-orbit are 
schematically drawn in Fig. 1. 

In the present framework, each valence-neutron orbit is introduced to have a definite 
value at the zero limit of centers of local Gaussians {{R^}) describing the spatial distribution 
of the orbit. The precise positions of {R^} are determined variationally before the angular- 
momentum projection. Since the values of {R^} are optimized to be finite, the orbits are 
not exactly the eigen-state of K^^ and are labeled as K^. For the vr-orbit with = 3/2~ 
(|3/2~)), the spatial part and the spin part of K are defined to be parallel (rYii|n |)), for 
which the spin-orbit interaction acts attractively. At the same time, |3/2~) is described as 
a linear combination of two orbits centered at right- and left-hand side of the system based 
on the MO picture. 

|3/2;) = -^{^{p^ + ipy)+a + -^{Px + iPy)-a}\n T). (5) 

{Px)±a = G±aez+bex ~ G±aez-hex) {Py)±a = G^±ae^+6e^ ~ G^±ae^-6ey- (6) 

Here, {px + ipy)±a denotes the p-orbit centered at ±a on the z-axis, and these variational 
parameters a and b are optimized by using the Cooling Method in antisymmetrized molecular 
dynamics (AMD) [|1^,|1^ for each basis state. Furthermore, the |l/2~) orbit where the spin- 
orbit interaction acts repulsively can also be defined by changing the spin-direction of |3/2~), 
where the spatial part and the spin part of K are anti-parallel (rYii|n |)): 

11/2;) = -^^^^^P^ + ^Py^^^ + -^^P^ + ^Py)-a}\^ i)- (^) 

The distribution of the a-orbit is just along the 3a axis, then it is introduced to have three 
nodes. |l/2~) is represented as a linear combination of three orbits with = l/2~, whose 
centers are +a, 0, and —a on the z-axis. 

|l/2;) = -^{iPz)+a - (p.)o + ip.)-a}\n T), (8) 



{Pz)±a — Gaez+bez " Gae^-hcz, {Pz)o — G^e^ — G-^g^. (9) 

These three orbits (|3/2~), |l/2~) and |l/2~)) are the basic building blocks for the molecular- 
orbital structure. Also, | — 3/2~), | — l/2~) and | — l/2~) orbits are introduced by taking 
the time reversal of |3/2~), |l/2~), and |l/2~) orbits, respectively. 

The Hamiltonian and the eflFective nucleon-nucleon interaction are the same as in Refs. 



Tq] , and parameters of Volkov No. 2 [^] for the central part and the G3RS spin-orbit term 



2T[| for the spin-orbit part are determined from the a + n and a + a scattering phase shifts. 



and the binding energy of deuteron is also reproduced with these parameters. 
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III. RESULTS AND DISCUSSIONS FOR BENDING STABILITY 



In the following part, we show the calculated results for the stability of the linear- 
chain state for various configurations. The isotopes and configurations which we take into 
account are ^^C, ^^0(3/2-)^ (two n's in the 7r-orbits), ^^0(1/2")^ (two n's in the a-orbits), 
^^C((3/2-)2(l/2-)2) (four n's in the 7r-orbits) and ^^C{{?>/2-f{l/2-f) (two n's in the tt- 
orbits and two n's in the a-orbits). As schematically shown in Fig. 2, two variational paths 
are introduced corresponding to the breathing-like (Fig. 2 (a)) and the bending-like (Fig. 2 
(b)) degrees of freedom. The parameters d and 9 stand for the a-a distance and the bending 
angle of the 3a-core, respectively. 

Firstly, we show the 0^ energy curves for the linear-chain structure against the breathing- 
path in Fig. 3. It is found that the energy pocket around = 3 fm becomes deeper 
as the increase of number of valence neutrons in the vr-orbit (-^^C ^"^0(3/2")^ 
^^C((3/2~)^(l/2~)^)). Using our framework, the binding energy of one o^-cluster is cal- 
culated to be 27.5 MeV, then the a+a+a +n-neutrons threshold energy is —82.5 MeV. The 
3Qf-system (^^C) has minimal energy around d = 3.5 fm, however, this is too shallow to con- 
clude the stability of the linear-chain state. On the contrary, in ^^C(3/2~)^, there appears 
evident minimal energy around d = 3 fm. The energy (- -82 MeV) is lower than by 
11 MeV and the energy pocket is much deeper. This energy corresponds to the excitation 
energy of 18 MeV from the ground state calculated with an equilateral-triangle configura- 
tion for the 3a-core, which is -101.2 MeV. The ^^C{{?>/2-f{l/2-f) configurational state is 
most stable among states studied and has an energy pocket of ^ —86 MeV, where the a-a 
distance is = 2.5 fm, shorter than those for -^^C and -^^0(3/2")^. Therefore, the vr-orbit 
is found to stabilize the linear-chain structure as the increase of valence neutrons (^^C 

Next, we discuss the case where the valence neutrons occupy the a-orbit. ^^C(l/2^)^ 
has an excitation energy higher by about 14 MeV in comparison with ^^C(3/2~)^. It is 
rather surprising that the difference is only 14 MeV in spite of the fact that 3/2~ has 
only one node and l/2~ has three nodes. This is because the a-orbit is along the a-a-a 
core: The higher nodal orbits along the symmetry axis become low-lying as a result of the 
clustering of the core. The a-orbit enhances the prolonged shape of the 3a-core, and the 
optimal rf-value is ^3.5 fm. However, the -^^0(1/2")^ case has no deep pocket enough to 
be stabilized against the breathing-like path. When two more valence neutrons occupy the 
TT-orbit, although this minimal energy is higher by 5 MeV than that of ■^^C((3/2~)^(l/2~)^), 
^^C((3/2-)^(l/2-)^) has the minimal 0+ energy of - -77 MeV. The calculated energy 
pocket is deep enough to guarantee the stability for the breathing-like path. Therefore, it is 
summarized that the linear-chain structure cannot be stabilized against the breathing-like 
mode with the a-orbit only in -^^C((l/2~)^), however, the neutrons in vr-orbits increase the 
binding energy and prevent the break-up of the system in -^^C((3/2~)^), -^^C((3/2~)^(l/2~)^), 
and ''Cii3/2-ni/2-f). 

Finally, the stability of these linear-chain states against the bending-like path is ex- 
amined. The 0+ energy curves of ^^C, ^'^C{{3/2-f), ^'^C{{l/2-f), i^C((3/2;)2(l/2;)2), 
and ■^^C((3/2~)^(l/2~)^) against the 0-value are shown in Fig. 4. Except for the case 
of ■^^C((3/2~)^(l/2~)^), the curvature of these states is rather monotonic and the energy 
minimum does not clearly appear. In ^^C, the orthogonality between the linear-chain con- 
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figuration and low-lying states with equilateral-triangle configuration of 3a is taken into 
account. This effect works but not sufficiently to push up the energy of the state with a 
finite bending angle. However, as clearly seen in Fig. 4, the ^^C((3/2~)^(l/2~)^) case shows 
a sharp increase of the 0"^ energy as the increase of the bending angle and is found to be 
stable against the bending-like path. This feature is much different from ^^C, ^^C((3/2~)^), 
^^C((l/2^)^), and ^^C((3/2~)^(l/2;^)^) cases. From the analysis above, the linear-chain 
configuration can be stabilized against the breathing-like path by neutrons in the vr-orbit 
(^4C((3/2-)^), ''C{(3/2-ni/2-n ''C{i3/2-ni/2-f)), but i«C((3/2;^)^(l/2-)2) is only 
the case which is stable also against the bending-like path. 

We further discuss the reason for the sharp increase of the 0"^ energy in 
^^C((3/2^)^(l/2~)^) against the bending-like path. The 0"*" energy increases by 15.7 MeV 
from 9 = 0^ to 6 = 30^, in which the kinetic energy part is 10.3 MeV. To understand the 
energy increase with the increase of bending angle 9 of this case, we calculate and compare 
the overlap between the wave functions with ^ = 0^ and 9 — 30^ for various configurations 
in Table I, where the results are shown. In ^^C, the wave functions with ^ = 0^ and 9 — 30^ 
have the squared overlap of 0.91, and ^^C((3/2~)^) has almost the same value. ^^C((3/2~)^) 
has the value of 0.85, smaller than ^^C{{3/2-f) by only 6 %, and i^C((3/2-)2(l/2-)2) has 
almost the same value as the ■^^C((l/2~)^) case. This result shows that the overlaps addi- 
tionally decreases a little for the cr-orbital neutrons, and also for the 7r-orbital neutrons as the 
increase of the valence neutrons. In spite of these, the overlap between ^ = 0^ and 9 — 30^ for 
i^C((3/2-)2(l/2-)2) case shows a significantly large decrease to 0.60. i^C((3/2-)2(l/2;)2) 
is only the configuration which shows drastic decrease of the overlap between 9 — 0^ and 
9 = 30^. 

As discussed in following part, it can be known that the drastic decrease as the increase 
of the bending angle is due to the increase of overlap between two neutrons in the 7r-orbit 
and two neutrons in the a- orbit. When there arises the overlap between them, the overlap 
component in the total wave function is diminished due to Pauli exclusion principle, that 
is, so-called Pauli blocking. Therefore, the physical state can be expressed by the modified 
wave function which is made by subtracting the overlap component from the original wave 
function. Since the energies of the tt- and the a- orbits discussed here are relatively low, the 
modified wave function involves larger components of higher excitation energy in comparison 
with the wave function at ^ = free from the Pauli blocking. As a result, the Pauli 
blocking due to the increase of overlap between valence four neutrons is considered to bring 
the increase of energy proportional to the decrease of squared overlap. This is a possible 
explanation for rapid increase of the energy against the bending angle. On this view we 
can estimate the increase of the energy as follows. Firstly, we compare the cases with Pauli 
blocking and without Pauli blocking. If there is no Pauli blocking, the squared overlap of 
i^C((3/2-)2(l/2-)2) between ^ = 0^ and ^ = 30^ is estimated as follows: 

PC^Cj^a)NoPauli = — ^ ^r)fl2n\ ^ ^ 0-85. (10) 

Here, P(^^C^), P(^^C^), and P{^^C) are squared overlaps between = 0^ and 6^ = 30^ 
of i^C((3/2-)2), ^^C{{l/2-y), and shown in Table I. In the actual squared overlap of 
^^C((3/2~)^(l/2~)^) shown in Table I, the Pauh blocking effect is automatically included, 
and the value is PC'^C^a)withPauii = 0.60. This difference corresponds to the component 
excited to higher shells because of the Pauh blocking effect, and this is 25 %. 
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PC^C^a)NoPauli " PC^Gj:(7)withPauli — 0.25 (11) 

Secondly, we estimate the increase of the kinetic energy when two valence neutrons in the 
TT- or the (j-orbits are excited to higher shells by 25 % due to the Pauli blocking effect. 

AE^2xhujx 0.25 lOMeV. (12) 

This value is consistent with the increase of the kinetic energy at 6^ = 30^ mentioned above. 
Therefore, we can conclude that this is one of the most important mechanisms which stabi- 
lizes the linear-chain state in ^^C. 

IV. SUMMARY AND DISCUSSIONS 

It is summarized that the linear-chain structure of ■^^C((3/2~)^(l/2~)^) with 3a core is 
the only case to have the simultaneous stabilities for the breathing-like break up path and for 
the bending-like path among ^^C, -^^C, and ^^C. Other configurations, such as ^^C((3/2~)^) 
and ■^^C((3/2~)^(l/2~)^) are stable against the breathing-like path but not stable against 
the bending-like path. A combination of the tt- and the a- orbits occupied by four neutrons 
plays doubly important roles to make a deep energy pocket for breathing-like path and to 
prevent the bending-like free motion of the system. The ^^C((3/2~)^(l/2~)^) configuration 
forms a rotational band with an energy slope of |y = 150 keV, and since the two-neutron 
separation-energy of ^^C is experimentally found to be 5.5 MeV, the band head energy is 
expected to corresponds to around 25 MeV in excitation. 

We have shown in ^^C that the Pauli-blocking effect among the valence neutrons play an 
important role for the stability of the linear-chain configuration. This effect is expected to be 
more important as the neutron number increases. We are interested in ^^C the appearance 
of similar linear-chain structure in low-lying region, where all of the three neutron orbits 
introduced here are occupied ((3/2-)2(l/2-)2(l/2;)2). 

Experimentally, the observation of the 7-ray from high-spin states is the most probable 
way to confirm the presence of rotational band with large moment of inertia. Although the 
present analysis has been restricted to the 0"^ state, we will investigate also high-spin states, 
where it is expected that a glue-like role of valence neutrons becomes more important to 
prevent the breathing-like break-up of the system. 

In the present analysis, parameter d has been fixed when the stability against the bending 
motion has been examined, and 6 has been fixed to zero when the bending motion has 
been examined. However, to assure theoretically the stability of these linear-chain states, 
we intend to investigate the case when two degrees of freedom {d and 6) simultaneously 
activate. This shall be examined by superposing states on the energy surface of d and 9 
based on generator coordinate method. 

The authors thank members of RI beam science laboratory in RIKEN for discussions and 
encouragements. One of the authors (N.I) thanks Prof. R. Lovas, Prof. W. von Oertzen, 
Prof. H. Horiuchi, and Dr. Y. Kanada-En'yo for fruitful discussions. 
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FIGURES 



FIG. 1. The schematic figure for the amplitude of the 7r-orbit (a) and the cr-orbit (b) on the 
X — z plane. The circles represent the a-clusters. 

FIG. 2. The schematic figure for the breathing (a) and the bending (b) motion of the lin- 
ear-chain state. The stability of the linear-chain state is examined for the a-a distance {d in (a)) 
and bending angle {9 in (b)). 

FIG. 3. The 0^ energy curves against the a-a distance (d) for ^'^C (solid curve), ■^^C((3/2~)^) 
(dashed curve), ^'^C{{l/2-f) (dotted curve), ^^C{{?>/2-f{l/2-f) (dash dotted curve), and 
■'^^C((3/2~)^(l/2~)^) (dash two-dotted curve). The coefficients of local Gaussians {gj} describ- 
ing the valence neutrons are treated as variational parameters to take into account deviations of 
the original MO. 

FIG. 4. The 0^ energy curves against the bending angle {9) for ^'^C (solid curve), ■'^^C((3/2~)^) 
(dashed curve), ^^C{{?>/2-f) (dotted curve), ^^C((3/2-)2(l/2-)2) (dash dotted curve), and 
^^C((3/2~)^(l/2~)^) (dash two-dotted curve). The coefficients for the linear combination of Gaus- 
sians describing MO are optimized at ^ = 0^ and fixed. The a-a distance [d) is fixed to 3 fm. 
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TABLES 



TABLE L The squared overlap between ^ = 0" and ^ = 30° for ^^c, ^'^C{{^/2-f), 
C((3/2-)2), i6c((3/2-)2(l/2-)2), i6c((3/2-)2(l/2-)2), 



configuration 


squared overlap between = {f and 9 = 30° 


12c 


0.906 


i4C((3/2-)2) 


0.905 


i4C((3/2-)2) 


0.804 


i6C((3/2-)2(l/2-)2) 
i6C((3/2-)2(l/2-)2) 


0.865 


0.602 
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Fig. 3 
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